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I. INTRODUCTION 

Doped manganites (-Ri_ x , A x )„ + iMn„03 n+ i (R=La, 
Pr, Nd, Sm ; A= Ca, Sr, Ba ; n = l,2,oo) have re- 
cently attracted considerable interests due to the colos- 
sal magnetoresistance (CMR) observed near the rfeE|ro- 
magnetic (spin-F-type) transition temperature TjjtJ It 
is now recognized that the most fundamental interac- 
tion in these materials is the double exchange inte*ac«j 
tion, which connects the transport and magnetism.LTlLHI 
Since the discovery of the CMRej, however, it has been 
pointed out that the double-exchang-p mechanism alone 
cannot explain not only the CMREil but also the sev- 
eral observed properties in this system. As the mech- 
anism which plays an essential role on CMR in coop- 
erate with double-exchange interaction, several candi- 
dates has been, suggested, for example, the JahitTeller 
(JT) polaronliirtia, thcj-charge inhomogeneityliJtfl, the 
percolative processeai3'E3, the phase segregation with re- 
spect to the, otbitalEj, and the orbital polarization and 
fluctuationE3~Ea, the last of which we describe in this 
manuscript. 

The well-known issue showing the importance of the 
orbital polarization is the layered (A-type) antiferromag- 
netism observed in the mothec compound of this 

system.EjED Kugel and KhomskiiEa treated this in a 
framework of the superexchange interaction and showed 
that the full consideration of the orbital degeneracy is in- 
dispensable to explain the spin A-type structure. There 
the orbital polarization is essential: Under the doubly de- 
generate orbitals, the on-site Coulombic repulsion differs 
depending on the configuration of the occupation, U for 
the two electrons occupying the same orbital, U' — J for 
occupying the different orbitals with the parallel spin, 
and U' + J for occupying the different orbitals with the 
anti-parallel spins, where U and U' are the intra- and 
inter-orbital Coulombic interactions, respectively, and J 
is the interorbital exchange interaction.^ In order to 
maximize the energy gain via the second-order pertur- 
bative processes, electrons form the staggered orbital oc- 
cupation (AE. orbital ordering) with the energy gain by 
i 2 J(U' — J).t3 In such an orbital ordering, there is a def- 
inite distinction between the occupied and unoccupied 
orbitals, which is the orbital polarization. The orbital 
polarization (or orbital ordering) is the important origin 
of the j4-type spin structure in the mother compound. 
In the viewpoint that the CMR with x ~ 0.1750 occurs 
in the lightly doped Mott insulator, the orbital polariza- 
tion is likely to survive and to play an important role on 
CMR. 



Another point is that the origin of Hund's coupling 
Jff is nothing but the on-site Coulomb interactions. It 
seems therefore rather artificial to take Jh — ► oo while 
the on-site repulsion is neglected^as in the framework of 
the double-exchange mechanism. 

In a form which include the mother compound, we 
studied the extended Hubbard-type model witb-,|tiie or- 
bital degeneracy for any doping concentration.c3E.il Cal- 
culated meanfield phase diagram well reproduced the 
global topology of the magnetic structure depending on 
the doping concentration x.OO With a large orbital po- 
larization we could predict the emergence of the ^4-type 
and the rod-type (C-type) AF in the moderately doped 
region, independently from the experiments discovering 
these phases in the CMR compound withpfinitjC x (which 
is larger than that for the CMR region)! 2 ^^ 2 ] It turned 
out that these phases couldn't be reproduced without a 
large orbital polarization, which is therefore essential in 
the doped region where these AF phases are observed. 
Because the orbital polarization increases as x decreases, 
this concludes that the large orbital polarization survives 
even in the CMR region with smaller x. Base on this 
result, we discuss the spin canting, the spin wave disper- 
sion, and spin wave stiffness from the standing point of 
the large orbital polarization. We also discuss the or- 
bital fluctuation which turned out to be important in 
the ferromagnetic metallic (FM) region where the CMR 
is observed. Spin wave softening near the zone boundary 
is also discussed in this context. 



II. MODEL AND FORMALISM 



We start with the Hamiltonian 



H = H 



K 



H 



Hund 



H orl site + Hs + -ffcl-ph , (1) 



where Hk is the kinetic energy of e g electrons, -ffHund 
is the Hund's coupling between e g and t2 g spins, and 
H on s ito represents the on-site Coulomb interactions be- 
tween e g electrons. ti g spins are treated as the local- 
ized spins with S = 3/2. The AF coupling between 
nearest neighboring ti g spins is introduced in H$ to re- 
produce the NaCLiype (G-type) AF spin ordering ob- 
served at x = 1.0.E3 Using an operator d\ which cre- 
ates an e g electron with spin a (= ] , [) in the orbital 7 
[= a(d x 2_ y 2) 7 b(d 3z 2_ r 2)] at site i, each term of Eq. (|l|) 
is given as 



H 



K 



E 



(2) 



1 



Hund 



and 



(ij) 



(3) 



(4) 



t77 in Hk is the electron transfer integral between Mar- 
ty u 

est neighboring sites and it depends on th£-pair of orbitals 
and the direction of the bond as follows E3 



/ _3 V3 
t 11 = tnl 4 4 



/ _3 _v/3 
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and 





i i+z — L I q _2 



(5) 



(6) 



(7) 



to is the electron transfer integral between d 3z 2_ r 2 or- 
bitals along the z direction. The spin operator for the e g 
electron is defined as S £g i — \^ di ia <Japd%yf) with the 

Pauli matrices a. St 2g i denotes the localized tig spin on 
the i site with S = 3/2. H on s ite consists of the intra- 
and the inter-orbital Coulomb .uiteractions and the inter- 
orbital exchange interactionpj'EJJ 

Hon site U ^ ^j'tT^J'tJ- 
31 

+U' n jaa n jb<y , 

jcrcr' 

+ J dj a(T dj brT ,dj aa idjba 

= -E(^+^) ! (8) 



where rij ia — dj dj ia and n J7 = X) n i70"i an d t ne 

cr 

isospin operator describing the orbital degrees of freedom 
is defined as 



(9) 



Coefficients of lip-spin and isospin operators, i.e., a and 
(3, are given asoEZI 



a = U-^>0, 



(10) 



(11) 



The magnitude of the meanfield solution of the isospin 
operator (T) gives the energy splitting between the occu- 
pied and unoccupied orbitals, namely the orbital polar- 
ization. Therefore, the minus sign in Eq. (^|) means that 
the on-site repulsion in this system induces not only the 
spin polarization but also the orbital polarization as the 
interplay with the orbital degree of freedom. By this or- 
bital polarization, the anisotropy of the e g orbitals is fully 
reflected to the transport and hence introduces the lower 
dimensionality even in the system with the isotropic crys- 
tal structure (n = oo ; 113-system). The parameters 
a,/3,to, used in the numerical calculation are chosen as 
to = 0.72 eV, U — 6.3 eV, and J-^-J.O eV, being relevant 
to the actual manganesapxides.c3o The electron-phonon 
interaction is given as,tj 



H, 



el— ph 



(12) 



where g is the coupling constant and r (ili) is the mag- 
nitude (direction) of the lattice distortion of the MnOv 
octahedra. Values of r and v are taken from the observed 
elongation as, r ~ 0.028 and v = ±(V3/2)x - (l/2li 
(staggered as dj, x 2_ r 2 / d^ y 2_ r2 ) for LaMnC>3 (n = oo),E3 
and r ~ 0.01, v//z (elongatipa along c-axis) in (Lax_ K , 
Sr x ) 3 Mn 2 07 (0.3 < x < 0.4) B 

In the path-integral representation, the grand partition 
function is represented as 



S = J Y\ D St2 g iDdi 1 <,Dd n<J eyL-piy- J dri(r)| 



(13) 



with 



L(t)=H (t) + di^a (r) (d T - n) d^ a (r) , (14) 



where r is the imaginary time introduced in the path- 
integral formalism, and c?i 7(T , di ia are the Grassmann 

variables corresponding to the operators d\ and di ia , 
respectively. By introducing two kinds of auxiliary fields 
corresponding to the following mean-field solutions, til 



0S = (S e 
0T = (?) , 



Jh 
25 



St 2g I ■ 



(15) 
(16) 



we obtain the effective action with respect to these aux- 
iliary fields and S t2g , after integrating over the fermion 
variables as,EJ 



and 



D{<p}e 



SeffM 



(17) 



2 



(18) 



(ij) 



(19) 



kk' ;nn' 577' ;a/3 

= (— iu n — fl) 5kk'\nn'-,*r~t'\ot0-\-Mkk'-,nn'\~i~t'-,ctl3 , (20) 



Mkk> 

= £fe 7 5kk'5nn'5 a (3 



(21) 



1 



ik-Rj —iu) n T 



, (22) 



for x — S,T. 

In the meanfield approximation, the free energy is 
given as 



Fj 



MF - 



-ksT ■ Seff[<p°] + fiN , 



(23) 



where ip c denotes the saddle point of <ps,T- We seek the 
saddle point within the several assumed ordering con- 
figurations, as following: We consider four kinds of the 
spin alignment in the cubic cell: spin F, A, C and G 
(NaCl-type). For spin A, we also consider the possibil- 
ity of the canting characterized by an angle n which is 
(it) for spin F (A). As for the double-layered com- 
pounds (n = 2), we consider an isolated double-layer, 
for which the Brillouin zone contains only two fc-points 
along c-axis, because the exchange interaction between 
two double-layers is reported to be less than 1 / 100 com- 
pared with the intra double-layer one.Lj As for the orbital 
degrees of freedom, we consider two sublattices /, and //, 
oigach of which the orbital is specified by the angle 8iji 



In the random-phase- approximation (RPA), we ex- 
pand S e ff[tp] with respect to the small fluctuation dips 
from its mean-field solution ip s for the spin degrees of 
freedom, 



0s = 0s + s 0i 



(25) 



Denoting the perpendicular (parallel) component to the 
mean- field as 7? (<?), 

5tps(k, uj n ) = a(k, w„) + n(k, uj n ) , (26) 

the deviation of the action can be written asEl 

5S eS = ^K n (q, fi) 7T {qs + q, H)-7T (-qs~q, -fi) 

+J2 K M n ) "(Qs+q, n)-{nX7f(-q s -q, -O)}. (27) 

q,Q, 

Because the spin wave is the Goldstone boson, the con- 
dition 1^(0,0) = ,-Kx(0,0) = 0, can be derived. 
Coefficient of the diagonalized quadratic form is ob- 



where we have introduced the momentum representation, tained as K 



TU) 



Kjr ± iK x , zero-point of which 



(ifm)(<f, f2 — —iuo) = 0) gives the dispersion relation of 
the excitation w = w(q). K a (q, Q) can be expanded as, 



K a (q, il) 



*A-in+ J2 c a ql 

a—x,y,z 

Bn 2 + J2 c a q 2 a 

a—x,y,z 



• ■ Spin F 
■ Spin AF 



(28) 



where a = 1 (— 1) corresponds to spin up (down), re- 
spectively. We evaluate only the static spin-wave stiff- 
ness C Q — C a (x) because the dynamic spin wave veloc- 
ity evaluated by using the above expression inherently 
gives a misleading estimation; For the half-filled insu- 
lator, x = 0, it does not reduce to the energy order 
as the superexchange interaction ~ t 2 /U, giving rather 
the order of perhaps due to the inherent fault of 
the RPA. For the metallic region, x ^ 0, we cannot 
reproduce the correct dispersion-relation, because the 
Landau-damping is not properly treated in our calcula- 
tion where the Brillouin zone is discretized and thus the 
gapless individual-excitation is not correctly evaluated. 
C Q = C a (x) roughly reflects the exchange-interaction de- 
pending on x. 



71,11 



•)+' 



7 1, II 



(24) 



We also consider four types of orbital-sublattice ordering, 
i.e., F-, A-, C-, G-type in the cubic cell. Henceforth, we 
often use a notation such as spin A, orbital G (#r,0jj) 
etc.. Denoting the wave vector of the spin (orbital) or- 
dering as qs (<Tt), the ground state energy is given as a 
function of the spin ordering (n, qs), the orbital ordering 
(@i,iii <7t), and the lattice distortion (<?, r, v). 
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III. ORBITAL POLARIZATION AND 
FLUCTUATION 
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FIG. 1. Mean field phase diagram as a function of the 

carrier concentration (a;) and the antiferromagnetic interac- 
tion between ia_Kpins (Js) for the cubic system (n = oo 
; 113-system)EfH Dotted line (J s =0.009) well reproduces 
the change of the spin structure experimentally observed. 



kinetic energy gain non-monotonically via the change in 
the density of states with the van-Hove singularity.Ea This 
can also be ^^igin of the instability toward the phase 
segregation,E31l3'EJ though it does not occur in our cal- 
culation. 

With no orbital polarization- such a non-monotinic be- 
havior cannot be reproduced pi as shown in Fig. |[ This 
is because the anisotropy of the degenerate orbitals are 
mixed to disappear with no polarization (in this case we 
cannot say which orbital is occupied because of the hy- 
bridization) . 



Spin F a =8.1 , 0=6.67 Spin A a =8.1 



=6.67 
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Fig. [l] shows the zero-temperature meanfield phase 
diagram of the cubic system (n = oo ; 113-system) in 
a plane of x and Js (AF superexchange interaction be- 
tween t2 g spins), with the. optimization of the orbital at 
each point on the plane.c§E-3 With Js being fixed to bfip. 
relevant value to the actual compounds, Js =0.009,E3'EZl 
we obtain the spin transition as A — ► F — > A — > C -^r-G 
with increasing x, being consistent with experimentsca. 
Non-monotonic phase boundaries are essential for these 
variety of the spin structures. 




FIG. 2. Mean field phase diagram with no orbital polar- 
ization. In this case the nonmonotonic behavior of the phase 
boundaries disappears. 



Dimensionality control by the orbital polarization is the 
origin of such a behavior: Orbital ordering changes 
from that maximizing the superexchange energy gain 
for smaller x (< 0.3) to that maximizing the double- 
exchange energy gftin for larger x, with the change in 
the dimensionality.^ This orbital transition varies the 



FIG. 3. The energy as a function of the orbital state charac- 
terized by 9 in the several value of x. (a) Spin F is assumed, 
(b) Spin A is assumed. In both cases, the orbital F"-type 
structure is assumed.E-1 



For the global topology A ^ F ^ A — > C G to be 
reproduced, it is therefore necessary that a large orbital 
polarization occurs even in the spin F phase where CMR 
is observed. As for the spin F (CMR) phase, however due 
to its isotropy, the question is how does it coexist with, 
the observed isotropic properties in CMR compounds J3'E3 
because such a polarization leads to the anisotropic car- 
rier hopping. The key for this question is the orbital 
fluctuation. 

Fig. [| shows the energy dependence ou the orbital 
configuration for spin F and A phasesEJ In spin F 
phase, there are many degenerate saddle points due to 
the isotropy, and the height of the barrier is an order 
smaller than that for the AF phase. 
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FIG. 4. Schematic picture of the orbital liquid state. 



This result implies that the orbital fluctuation becomes 



4 



important in the spin F phase comparing with the AF 
phases. Reentrant of the spin A with increasing x, seen in 
Fig. [l] implies that the d x 2_ y 2 orbital ordering is inherent 
property of the double-exchange interaction. Therefore, 
in the extent beyond the mean field theory, it is likely that 
the degenerate saddle points, d x 2_ y 2 1 d y 2_ z 2, and d z 2_ x 2 
resonate to recover the isotropy of the spin F metallic 
phase though the large orbital polarization still survives, 
forming the orbital liquid state, as shown in Fig. 



IV. SPIN CANTING AND ORBITAL ORDERING 

Fig. U shows the phase diagram of the layered com- 
pounds (n = 2 ; 327-system)£3 In this system, the 
anisotropic crystal structure also controls the dimension- 
ality: restricted hopping along c-axis brings about the 
d x 2_ y 2 -orbital ordering in the metallic region even for 
the isotropic spin G- and F (x > 0.2) alignment. Es- 
pecially the planer spin F phase seen for x > 0.2 js-iQS-i 
sential for the spin canting observed in this systemE§CJ 
with 0.4 < x exp , < 0.48, as below. The global topol- 
ogy of the phase diagram in this system is reproduced as 
A — > F — > A — > G with increasing xjt3 as shown in Fig. 

where there are two phase boundaries between the 
spin A and F phases; One is with small x (x < 0.1, left 
boundary), and the other is with finite x [x > 0.1, right 
boundary). Under the competition between the super- 
and double-exchange interaction, the saddle point of the 
canting angle rj is given as,L3 



V 



t z x 
4Js 



(29) 



where the t z denotes the inter-layer hopping integral. 
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i i i i F^" 

0.1 x yo.2 \ 0.3 

/ p-spin F 

orbital transition 

FIG. 5. Phase diagram for the layered compounds (n = 2 

; 327-system)£3 



The right-hand-side of the above equation should be 
smaller than the unity for the occurrence of the spin cant- 
ing. For the left boundary, this condition can be satisfied 
for any magnitude of t z because the small x always makes 
the right-hand-side of Eq. (|2^) to be small. For the right 
boundary, however, this can be satisfied only when the 
orbital is planer, i.e., nearly d x 2_ y 2 with small t z ; the 
characteristic energy scale of the hopping integral is an 
order greater than that for Js, therefore if the orbital is 
spherical, t z /Js becomes much more than the unity and 
hence the right-hand-side of Eq. fl29| ) because in this case 
the finite x does not make it small any more. Therefore 
it is concluded that the planer orbital is indispensable for 
the canting observed on the right phase boundary with 
finite x (metallic canting). 

Experimentally, this metallic canting is comraxml* 
found in the double-layered compound (n = 2).tffl 
This is because the layered structure stabilizes the planer 
(d x 2_ y 2) orbital in the metallic region. In the 113-system, 
on the other hands, no spin canting on the right bound- 
ary is reported,Ej which may be accounted by its isotropy 
leading to no such stabilization. This isotropy of the 113- 
compounds may allow only two possibilities for its orbital 
state; one is the orbital liquid state resonatingj-among the 
planer orbitals, d x 2_ y 2, d y 2_ z 2, and d z 2_ x 2fn and the 
other is the quasi-spherical orbital, which is obtained as 
the saddle point within thc-extent of the mcanfield the- 
ory, as shown in Fig. |l].oE-3 Kawano et al. observed 
the metallic canting in 113-system, Ndi^Sriy^MnOs, 
with slight anisotropy of the lattice structureCJ in the 
temperature-driven transition between the spin F (high- 
temperature phase) and the spin A (low-temperature 
phase). This supports the former possibility of the or- 
bital, i.e., the orbital liquid state; If the orbital is quasi- 
spherical in perfectly cubic system, taking the latter pos- 
sibility, such a slight lattice anisotropy leads to only a 
slight distortion of the spherical orbital which remains 
the right-hand-side of Eq. (|2^) still larger than the unity 
and hence no canting is expected. On the other hand, 
taking the former possibility, such a slight anisotropy is 
enough to stabilize d x 2_ y 2 immediately and hence the 
metallic canting can be explained. 

Another important feature as for the metallic canting 
is the stability of the spin A phase against the canting. 
When the x holes are introduced, the kinetic energy gain 
AEkin (£) via Uie bonding and anti-bonding splitting A = 



t z cos 4 — iz£.u is given as 



AEkin (£,) ' 



-tl-N F -e (for£<? c 
-t z ■ x ■ £ (for £ > £ c ) 



N F t z 



(30) 



with simplifications of a perfect spin polarization and 
the constant density of states. The competition between 
this kinetic energy gain and the energy cost of the ex- 
change interaction, Jscosry = Js (2£ 2 — l), is the origin 
of the spin canting. The lower line of Eq. ( pSO] ) is ob- 
tained hv de Gennes, and if this holds the canting always 
occursi!3 The new aspect here is that AEki n (0 oc £ 2 
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when the splitting A = t z £ is smaller than the Fermi en- 
ergy e F = x/N F and both the bonding and anti-bonding 
bands are occupied. Therefore the spin A structure 
(£ = 0) V = 7r ) is a t least locally stable when 2 Js > t 2 z Np. 
This condition can be satisfied when the orbital is almost 
d x 2_ y 2 and t z is much reduced from to- By minimiz- 
ing the total energy AE(£) = A£ Wn (0 + &Eex(£), it is 
found that the spin canting can occur only when £ c < 1; 
When £ c > 1 (x > t z Np), only the upper line of Eq. 
© is relevant and AE = (2J S - t 2 z N F ) ■ £ 2 . Therefore 
£ jumps from 1 (spin F) to (spin A) as Js increases 
across t 2 z N F /2. 



(spin F) ^=1 



z F 



(spin A) £=0 



- xt 2 
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FIGi-6. the optimized £ as a function of Js for the case 



When £ c < 1 (the spin canting can occur), the opti- 
mized £ as a function of J5 is given in Fig. 0. As Js 
increases, the spin structure changes as spin F (Js < 
t z x/A) — > spin canting {t z x/4 < J s < t 2 N F /A) — > 
spin canting with fixed canting angle (t 2 z N F /A < J? < 
t 2 z N F /2) -> spin A (t 2 z N F /2 < J s ). Note that the cant- 
ing angle continuously evolves from spin .F, but jumps at 
the transition to the spin A. This seems to be consistent 
with experiments^ where the canting angle larger than 
63 deg. is not observed. 



V. SPIN DYNAMICS AND ORBITAL 



By fitting K a (q, 0) as a function of q, in Eq. (f28|), we 
can evaluate the static stiffness of the spin wave excita- 
tion C a due to the e g orbital contribution. As the orbital 
configurations to be assumed, we take the saddje-point 
solution obtained in the meanfield theory as,r 3 



x = 0.0 Spin A Orbital C 

x = 0.1 Spin F Orbital C 

x = 0.2 - 0.4 Spin A Orbital F 

x = 0.5 - 0.9 Spin C Orbital F 



(60, -60) 
(80, -80) 
(0,0) 

(180,180). 



As for x — 0, we further introduced the JT effect[EI bj4 
putting the observed distortion of the MnOg octahedral^ 
Fig. m shows the g-dependence of — K\ (q x ,0) for the 



spin A configuration with d x i_ y i orbital ordering (Minus 
sign of comes from the negative B in Eq. ( |28|) to 
correspond the positive sign of the plot with the stability 
of the saddle point). 



0.08 



(100)-direction 



3 



0.04 



o 



-0.04 




% spin A, orbital F(0,0) 

FIG. 7. wave vector dependence of —Ky (q x ,Q) calculated 
for the spin A, d a .i_ y i orbital ordering, as an example. Minus 
sign of Ki comes from the negative B in Eq. ( j28|) to corre- 
spond the positive sign of the plot with the stability of the 
saddle point. 

We have chosen this structure, as an example, because 
the double exchange interaction is most-effective in this 
ordering as the meanfield theory shows ,E-3 and hence the 
crossover from super- to double-exchange manifests it- 
self most remarkably. The enhancement of the stiffness 
with increasing x can be reproduced. This is due to the 
crossover from the super- to the double-exchange inter- 
action as x increases, which is well evaluated in our for- 
malism in the unified way, as the inter- and intra-band 
transitions, respectively. Plots are well fitted in the whole 
Brillouin zone for the ferromagnetic-bond direction by 



K l (<f>0) oc (1 - cosq a ) , 



(31) 



not only in this case but also for all the other order- 
ing shown in the above table. This implies that only 
the nearest-neighbor interactions are important in the 
spin wave excitation. This issue is important because 
there is no guarantee that the exchange interaction can 
be represented by the nearest neighbor Heisenberg model 
at finite doping, and because the softening near the 
zone boundary—has been observed in some materials 
experimentallyo uj Our .result here is in sharp contrast 
to the first principle studyt£l which attributes the origin of 
such a softening to the longer-range interactions than the 
nearest- neighbor interactions. Negative stiffness C a < 
seen for x — 0, (lOO)-direction, corresponds to the insta- 
bility of the spin structure, which can be explained as fol- 
lows. Around x — the spin structure is dominated by the 
superexchange interaction where the energy gain for spin- 

F (AF) bond is t 2 _ u /p (t 2 _ /a)B where t _ (t _ u ) 
are the transfer integral between the nearest-neighboring 
occupied/occupied (occupied/unoccupied) orbitals. Or- 
bital F (0,0) leads to t x l > t*™ u and thus the intra- 
plane bonds favor spin-^J 1 for our choice of the param- 
eters a rs /3. This destabilizes the spin A structure in 
(lOO)-direction. As the doping x increases, the double- 
exchange interaction, becomes more and more important. 
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This stabilizes the ferromagnetic bond within the plane, 
and C x ( y \ becomes positive. 
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FIG. 8. Doping-dependence of the spin stiffness. The or- 
bital and the spin structure are optimized at each point. The 
enhancement of the spin-stiffness and thp-icross-over of the 
dimensionality are seen with increasing x.E3 



Fig. |g shows the static spin wave stiffness as a function 
of the doping concentration x, including the contribution 
from ti g (Js)- We could reproduce the qualitative fea- 
ture of the dimensional crossover and the enhancements 
of the stiffness in terms of the crossover from the super- 
ior smaller a;) to the double-exchange interactions (for 
larger x) accompanied with the change in the orbital or- 
dering. As x increases, the spin structure changes from 
spin A insulator at x = into the nearly isotropic spin F 
metal, to the spin A metal with two-dimensional d x 2_ y 2- 
orbital .alignment, and to the spin C metal with d 32 2_ r 2- 
orbital.E^IEJ Accordingly, the in-plane stiffness shows an 
increase, moderately at the beginning and then rapidly in 
the region of the spin A-metal. This reflects the fact that 
the double-exchange interaction is the most effective and 
prefers the d x 2_ J/ 2-orbital, i.e., the double-exchange inter- 
action is basically two-dimensional with the e ff -orbitals. 
In the spin C metal for x > 0.4, one-dimensional orbital 
along (OOl)-direction gives rise to a steep increase of the 
stiffness in this direction. 

The observed anisotropy of the spin stiffness is de- 
termined by the long range ordering of the orbitals. 
Fig. U also represents the cross-over of the dimcEt. 
sionality which we proposed in the previous report £a 
The stiffness changes from the nearly isotropic one 
in the spin F state to the considerably strong two- 
dimensional one for spin A metal, and to a quasi- 
one-dimensional one for spin C, reflecting the ptbital 
transition with increasing x. Yoshizawa et aLE3 ob- 
served such two-dimensional anisotropy of the stiffness 
for Ndo.45Sro.55Mn03, being consistent with our re- 
sult. Quasi-one-dimensiopaLanisotropy is predicted for 
Ndi-^Sr^MnOs (x > 0.6)BB 



The in-plane spin stiffness ^totals-total m Fig- || could be 
compared-with the experiments. In Lai-^Sr^MnOs, En- 
doh et alz3 observed the plateau of the velocity v x in the 
orbital-ordered insulating state up to x ~ 0.12 and then 
the velocity increases in the spin F metallic phase. Com- 
paring this with the calculation above, it seems that the 
moderate increase up to x ~ 0.15 in Fig. |^ corresponds to 
the plateau, while the rapid increase foip-s; > 0.15 to the 
increasing velocity observed by EndohO Then orbital- 
ordered spin F metallic state in Fig. || corresponds to 
the insulating spin F phase in experiments. Both the 
spin F- and A-metal in experiments, on the other hand, 
seems to corresponds to the spin j4-metal with d x i_ y i 
orbital ordering in the calculation. This fits well orbital 
liquid picture by Ishihara et alE3; In the perfectly cu- 
bic system the orbital state in spin F metal is described 
as the resonance among d x 2_ y 2, d y 2_ z 2, and d z 2_ x 2. In 
the actuaL-CMR compound, however, the slight lattice 
distortionclJEJ may breaks the cubic symmetry to sta- 
bilize d x 2_ y 2, though it is still accompanied with large 
fluctuation around it. 

Now we turn to the absolute value of the stiffness 
in the spin F metallic phase. Taking the reported lat- 
tice constant and the magnitude of spin moment as, 
S* = 3/2 + 1/2(1 — x), the experimental values of 
the static spin stiffness, ^tal^otal' are H-61 meV 
La .7Sro.3Mn0 3 El and 10.24 meV for Ndo.7Sro.3Mn.O3c 
respectively. These are in quite well coincidence with 
"^totai^totai =10-53 meV, estimated by RPA here with 
x = 0.3, d x 2_ y 2 -orbital ordering. A simple tight-binding 
estimation of the static spin stiffness, 



D = 



S* d 



% 2 ) £ Uj 



(okU- CT |o> 



(32) 



with d a; 2_ y 2-orbital also gives the similar value, ~ 10 
meV (with to = 0.72 eV, x = 0.3), where the strong 
Coulombic interactions are reflected as the full orbital 
polarization d x 2_ y 2 (superexchange interactions are ig- 
nored). This agreement can be understood in terms of 
the above orbital liquid picture as follows: While the 
large orbital fluctuation around d x 2_ y 2 may cause the 
several anomalous behaviors in the transport properties, 
it is not reflected to the stiffness constant because the 
correction due to such a fluctuation has the wave vector 
dependence as ~ (1 — cosq a )i!3 as described in the next 
paragraph, doing little around q = and hence the stiff- 
ness constant. Therefore the d x 2 _ y 2 -orbital ordering can 
give a good estimation of the stiffness constant of spin F 
metallic phase with a large orbital fluctuation. 

The softening observed near the zone boundary of 
the spin wave excitation pan be understood in terms 
of the orbital fluctuationE3 When the normal vector 
of the resonating planer orbitals, d x 2_ v 2, d y 2_ z 2, and 
d z 2_ x 2, points along some bond direction, the ferromag- 
netic double-exchange interaction disappears along this 
bond resulting, instead, the AF interaction due to tig 
orbitals. Such an interaction between the orbital fluctu- 
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no (canceled out) coVeet^fo7(i being'consis- X>gKI 2 )- Then the polaronic effect is to replace [7 



ation and spin degrees of freedom leads to, in the lowest 
order, the self-energy correction with the fc-dependence 
as ~ (1 — cosg Q ) 2 for (0,0, £)- and (0, ^^-direction but 
no (canceled out) correction 
tent with the experimentsE3. 

The important implication concluded from the agree- 
ment between the experimental and RPA-estimated value 
of the stiffaess constant is little influence of the JT 
polar on Jiirtfl at least on the spin dynamics. JT po- 
laron should reduce the double-exchange interaction in 
the doped region via a bandwidth reduction. To describe 
this polaronic effect, we introduce here a generic model; 
Assume that the orbital configuration is relaxed to its 
stable one when the electron is occupying the site i. We 
express the polaronic degrees of freedom by the bosons. 
Now the electron operators <v , d have no orbital index, 
because of the sufficient orbital polarization, 

F = £ 4r«*i* + £ £ 9q0> q + bl q )d\ a d l<7 

ij,cr i.a q 

n^mx . (33) 



U case, the integral is determined by the small r re- 



gion where (Xj (t)X 3 (t)X] (0)X t {0) 



-At 



(A 



by U + A in the expression for J, which is a minor 
correction when U >> AEa, being in sharp contrast to 
the double-exchange interaction discussed above. Pola- 
ronic effect should therefore correct the RPA-estimation 
of the stiffness-enhancement as x increases to be smaller. 
Agreement between the observed and estimated stiffness 
in the doped region implies therefore that the spin dy- 
namics is not so affected by the JT polaron. This is also 
pointed out by Quijada et alEB 

Because the estimation is made under the assumption 
that the orbital is almost fully polarized to d x 2_ y 2, the 
agreement also suggests the large orbital polarization. 
With the absence of the orbital polarization, the stiff- 
ness enhancement should scale to electron density (1 — a;) 
rather than the hole x. The observed stiffness enhance- 
ment with increasing x even in the metallic region there- 
fore also supports the large orbital polarization. 



This is the usual polaron Hamiltonian, and the following 
unitary transformation U eliminates the coupling terms 
between electrons and bosons, 



U = exp 



EE 



fct 



(34) 



In 

ad 



terms of this U, the Hamiltonian H is transformed 

H = U^HU 

= E UjXtXjAvdjv + £ iu q b% 



^2 A ■ n itT + U ^2 n i] n ii 



(35) 



where X % = fflq>E,e*«-*«(g,/w,)(&, - &!,)], and A - 
53 g 5q/ w g i s the relaxation energy. We now derive the 
exchange interaction between spins in terms of the per- 
turbative expansion in ty . The double-exchange interac- 
tion is the first order in , and is reduced by the factor 
of < XjX j >= exp[-£ g K| 2 /2] („, = (g q /u q )(e^ - 

e lq ' Rj )) 1 which is exponentially small when g q /to q is large. 
This factor is nothing but the bandwidth reduction fac- 
tor due to the polaronic effect. On the other hand, for 
x = 0, the superexchange interaction under the coupling 
with the polaron is given by, 

J = 4|ty \ 2 J Q drGl (r)(xj(r) X d (t)x] (0)^(0)) , (36) 

where Gq(t) = e~ Ur ' 2 is the Green's function for local- 
ized electrons. Because we are interested in the large 



VI. CONCLUSIONS 

We discussed the zero-temperature phase diagram and 
the spin dynamics of the CMR compounds based on the 
model with a large orbital polarization. The topology 
of the magnetic transition depending on the doping con- 
centration cannot be reproduced without a large orbital 
polarization. This is because the double-exchange in- 
teraction is the most effective and prefers the d x 2_ y 2- 
orbital, i.e., the double-exchange interaction is basically 
two-dimensional with the orbital polarization. As for the 
ferromagnetic metallic phase the large orbital polariza- 
tion recovers the isotropy of the transport by forming 
a liquid state, i.e., the resonance among d x 2_ y 2, d y 2_ z 2 
and d z 2_ x 2. Spin A phase seen in the moderately doped 
region has a stability against the canting with infinites- 
imal angle deviation from n, being in sharp contrast to 
the spin A insulator. Though it cannot be infinitesimal, 
finite canting angle between (spin F) and ir (spin A) 
can realize only if the orbital is planer both in spin F 
and A. The observed metallic canting in f f 3-compounds 
is therefore an evidence that the ferromagnetic metallic 
phase consists of such a planer orbital. The dispersion of 
the spin wave excitation evaluated in the RPA is well fit- 
ted by the cosine curve. This implies that the excitation 
is almost dominated by the nearest-neighbor exchange 
interaction even in the double-exchange regime, being in 
conflict with the first principle result. Estimated stiffness 
constant shows good agreement with the observed values 
for the metallic region. This strongly implies the absence 
of the JT-polaronic influence on the spin dynamics in the 
doped region. Based on the above orbital liquid picture, 
we could explain the spin wave softening near the zone 
boundary, its anisotropy, and no influence due to the or- 
bital fluctuation on the stiffness constant. 
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